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ABSTRACT 



A method for identifying and - after identification - 
reducing linear time invariant systems on the basis of input - 
output records of the system is reviewed, applied to physical 
devices and extended to handle the case where noise is present. 
The method is implemented using an analog to digital converter 
in order to have the input and output records in a digitized 
form and thus to be able to use a digital computer program 
composed of a numerical integration subroutine and another 
subroutine to solve overdetermined sets of linear algebraic 
equations. Several practical examples are presented to 
demonstrate the accuracy and present capabilities of the 
procedure . 
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I. 



INTRODUCTION 



A. THE IDENTIFICATION PROBLEM 

All the modern techniques of analysis and design of 
any kind of system are based on the previous knowledge of 
a mathematical model of the system to be studied. This 
principle is also applicable to control systems. The form 
of the model depends on the methods of analysis and design 
to be employed as well as on the physical characteristics 
of the system and since most of the theory on control 
systems is based either on the state space or transform 
representation of systems, the majority of the mathematical 
models will be a set of state equations or a transfer 
function. 

After deciding what basic form the mathematical model 
should take, the problem of determining the numerical values 
of the parameters arises. These values can often be deter- 
mined from the laws of physics and the data supplied by the 
manufacturer or obtained through testing. This is not always 
the case, however; sometimes the laws of physics become 
mathematically intractable or are very difficult to apply and 
quite often the values of certain key parameters are not 
available or have changed with time. It is in these cases 
when the identification problem arises. This problem consists 
in determining the system model given the input, output and 
initial conditions. 
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For the purpose of this work, the identification problem 
can be stated as follows: 

Given a system considered a black box, obtain a record 
of the input and output of the system over some finite 
period of time - 400 msecs, in this particular case - and 
from it find a mathematical model - the transfer function 
relating output and input in this case - and the values of 
the parameters in such a way that the model will accurately 
describe the system. 

It should be kept in mind, however, that the problem of 
identifying a system solely on the basis of input and output 
quantities is very rare in engineering, because in most of the 
cases the engineer has some idea whether the system is linear 
or non-linear, time varying or time invariant, what is the 
order of the system itself and perhaps what are the time 
constants involved and in the presence of noise the statis- 
tical characteristics of it. For these reasons father than 
having a BLACK BOX, what the engineers have is what could be 
called - a GREY BOX -, and as will be shown later, the 
identification techniques depend very heavily on knowledge 
of the characteristics mentioned above and in this particular 
case THE ORDER OF THE SYSTEM. 

Although the control systems literature on system identi- 
fication techniques is quite vast, there are no known identi- 
fication schemes capable of handling all identification 
problems effectively. Some of the proposed identification 
methods are mentioned in the following: 
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Nieman, Fisher and Seborg [1] have pusblished a paper 
summarising most of the common approaches to the problem of 
identifying lumped parameter systems. Eykhoff [2] has a good 
discussion of the industrial application of several methods. 

Since the response to an impulse gives directly the transfer 
function of linear systems, obtaining this response is a 
common scheme used in linear identification. Mishkin and 
Haddad [3] have developed a technique for finding the impulse 
response based on samples of the system output and its deri- 
vatives. Levin [4] and Kerr and Surber [5] have developed a 
scheme for estimation of the impulse response from the 
samples of a noisy input and output. Turin [6] and Lichtenberger 
[7, 8] have used a matched filter to obtain an identification. 
Hill and McMurtry [9] have suggested the use of white noise, 
a binary test signal and crosscorrelation techniques to 
attack the problem. Cooper and Lindenlaub [10] have investi- 
gated the use of matched filter and crosscorrelation. 

Another common approach to the problem is to determine 
the coefficients of the differential or difference equation 
which describes the system. Kumar and Sridhar [11] have 
employed the method of quasilinearization with some success. 
Nagumo and Noda [12] have developed a learning approach to 
the problem. Bass [13] uses modulating functions and the 
method works well in the presence of noise. Astron and 
Bohlin [14] use a statistical optimal method of determining 
differential equation parameters known as the "maximum 
likelihood method". Peterka and Smuk [15, 16] have developed 
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a similar method which is not optimal but is considerably 
simpler computationally and it is known as the "prior 
knowledge fitting method". Ho and Kalman [17] have proposed 
an algorithm for determining state variable models of sampled 
data systems which works well in the presence of noise, and 
has been extended to continuous systems by Eldem [18] . 

Identification methods vary widely with respect to how 
much must be known about the system before the method can be 
applied. In general, the less known about a given system, 
the more difficult it is to find a method capable of accom- 
plishing the identification. 

B. OBJECTIVES OF THIS PAPER 

This paper will present a research, handling actual 
hardware, of an identification technique originally suggested 
by Diamesis [19] and studied using simulation techniques by 
W. R. Hansell [20] ; also an attempt at using the technique 
in the presence of noise is made. The method requires a 
knowledge of the system input and output over some finite 
time interval and a rough estimate of the order of the system 
is desired. Although, theoretically the system input used 
need not be restricted to a class of testing functions, that 
is, it can be completely arbitrary, experimental results 
showed that inputs such as step, ramp or in general any 
periodic wave, sine, square, etc., should be avoided. 

Unlike some identification schemes which require the 
calculation of derivatives of the input and output, the 
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technique to be presented requires only integrals of the 
input and output, and the advantages of numerical integration 
over differentiation are well known. 

In Chapter II , the theoretical development of the 
identification technique is given. A method for identifying 
the initial conditions of the unknown system is also presented, 
but in this case it is necessary to know the value of the 
derivatives of the input at time zero. Since the initial 
conditions do not modify in any sense the mathematical model, 
their knowledge is not important. 

Chapter III presents the method used to check the 
developed theory in Chapter II. Several examples are given 
to demonstrate the capabilities of the method. 

In Chapter IV, recommendation and conclusions are given 
and following this chapter a complete listing of all computer 
programs are presented. 
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II. IDENTIFICATION BY MULTIPLE INTEGRATIONS 



A. GENERAL APPROACH 

The development that follows is similar to the ones 
given by Diamesis [19] and W. R. Hansell [20] . The develop- 
ment given by Diamesis is restricted to the case where all 
initial conditions are zero and the parameters of the model 
are the unknowns of a uniquely determined set of linear 
algebraic equations. Hansell extended the development to 
the case where the initial conditions need not be zero and 
he makes use of an overdetermined set of linear algebriac 
equations with the model parameters and some coefficients 
which have information about the initial conditions , as 
unknowns . 

The overdetermined set of equations are solved using 
the method of least squares. The method is extended to 
handle noisy inputs and outputs. 

Any single input, single output, lumped parameter, linear 
time invariant system can be described by a linear ordinary 
differential equation with constant coefficients. The 
general form of this equation is given in equation (1) . A 
set of initial conditions is also necessary to completely 
define the system. 

y n) (t) + a n _ 1 y n_1) (t) + + a Q y(t) = ' 

( 1 ) 

m) v 

= b^U (t) + + b Q U(t) where, in a general case 

n>m and 
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y^ (t) and U-^ (t) represent the i tb and j tb derivatives 
of the output and input, respectively. 

The necessary initial conditions are: 

The value at time zero, of the output and its n-1 first 
derivatives and the value at time zero, of the input and its 
m-1 first derivatives. 

The notation to be used is: 



y (o) , y^ (o) y 11 2 ^ (o) , y 11 ^ (o) to denote 

the, n-1, I.C. of the output and 

u(o), u^ (o) u ra 2 (o) , u m 1 (o) to denote the. 



m-1, I.C. of the input. 

The identification problem to be treated here consists of 
determining the coefficients a^ and b^ of equation (1) , on 
the basis of input and output records taken over some arbi- 
trary time interval, for every pair of values chosen for n 
and m. The initial conditions will be assumed to be unknown 
and also it is assumed that the noise present is zero mean 
and that the standard deviation is small compared with the 
R.M.S. of the input and output, for the technique to work 
properly. 

Taking the Laplace transform of (1) , yields equation (2) 

s n Y(s) + a 1 s n-1 Y(s) + + a Y(s) = 

n-1 ' o 

b m smU(s) + + b Q U(s) + d n _ lS n 1 + ... + d Q (2) 

where the d^_ coefficients account for the contribution of 
the initial conditions. 

n + 2 

Dividing equation (2) by s n results in equation (3) 
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Y(s) 



+ a 



Y(s) 

n-1 2 

s 



+ a *<£> = b -2i§> 
o _n+l m n-m+1 



L U(s) , 1 

••• + b o “H+T + d n-l“? + 



s 
+ d 



1 

° -n+r 



(3) 



•Taking now the inverse transform of equation (3), gives 
equation (4). 




And since the system is assumed to be time invariant, 
nothing has been lost by setting the lower limit on the 
integrals equal to zero. 

Rearranging terms in equation (4), equation (5) is 
obtained. 
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y (t) dt n+1 + . . . 





Since records of the input and output are assumed to be 
known, a linear algebraic equation with the system parameters 
a^, bj , and the terms that contain information about the 
initial conditions as unknowns, can be formulated by performing 
the multiple integrations, from zero to some time t^, or from 
any time t Q to a final time t^, after assuming some values 
for n, and m. A set of 2n+m+l equations can be obtained by 
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letting t^ take on 2n+m+l different values. Assuming that 
the equations are linearly independent, it will now be 
possible to solve for the 2n+m+l unknowns; n-Hrt+1 , a^, b^ 
coefficients and the n d, terms which contain information 
about the initial conditions. It will be shown later how 
those coefficients can be used to determine the initial 
conditions. 

Consider now the case when noise is present at the input, 
at the output, or at both. 



Under these conditions equation (5) becomes: 
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y(t)+v(t) dt 




where v(t) and r(t) represent the noise at the output and 
input respectively. 

If the coefficients in equation (6) , that is, the 
multiple integrals of the input, output, and time, are going 
to have the same values as the corresponding ones in equation 
(5), in order to have the same solution for the a^, b ^ , and 
d^,, in other words, the same identification, the presence of 
noise should not affect the value of the integrals. 

Consider a general term: 




Equation (7) can be written: 



(7) 




The first integral of equation (8). 





y (t)dt p 




v(t)dt p 



( 8 ) 



(9) 



Jo J° 

is the true value of the coefficient of a .. 

n-p+1 

(5) and 



in equation 
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r t 



( p ) ••• 



v (t) dt p 



( 10 ) 



o I o 

the undesired contribution of the noise. Thus, what has to 
be done is to find when the value given by equation (10) is 
negligible with respect to the value of equation (9) . 

Assuming zero mean noise, the expected value of 
equation (10) is: 





••• (p) ... 




v (t)dt P 



••• (p) • • • 




H 



dt P = 0 



( 11 ) 



o jo 

Since zero mean noise was assumed and the expected value 
of the integral is equal to the integral of the mean value. 

But the fact that the mean value of equation (10) is zero is 
not enough to say that the influence of the value of equation 
(10) in (7) is negligible. Besides the mean of equation (10) , 
its standard deviation or R.M.S. value should be found in 
order to be able to compare it with equation (9) , and if the 
R.M.S. value of equation (7) is very small compared with the 
value of equation (9), it can be said that the influence of 
the noise will not significantly affect the identification. 

The variance of equation (10) considering, p = 2, for 
simplicity, is given by equation (12) 
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' \ 



= E 





v (y)dyd£ 



( 12 ) 



I J o Jo ' 

where y and £ are dummy variables. Equation (12) can be 
written 




can be written as R (t) , where t =y^-Y. . 

v ’ 2 1 1 



The R.M.S. value of equation (10) can now be expressed 
as follows: 



a 



z 





Rv(T)dYidY 2 d?id5 2 



1 

7 



< 





/ 



R v (o)dyd£ 



R v (x)< R v (o) 



(13a) 



o 



o 
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or going back to the initial notation 



a 



z 




^/R v (o)dt 2 



(14) 



If now (14) is compared with (9), with p = 2, it is observed 
that both integrals have the same form. Therefore, if 

y(t) >> R v (°) = R.M.S. value of v(t), the noise is not 
going to have much influence in the identification. The 
same conclusion can be reached for any value of p. 

As was said before, the set of 2n+m+l linearly indepen- 
dent algebraic equations obtained from equation (5) letting 
tf take 2n+m+l different values, would solve the problem of 

identifying the parameters a., b., and d, , but since for 

1 D K 

some kind of inputs and some t^, the 2n+m+l equations may 
not be linearly independent, thus instead of forming just 
2n+m+l equations, a set is formed which has a number of 
equations substantially bigger than 2n+m+l and this over- 
determined set is solved using the method of least squares. 

This technique also allows some flexibility in the 
selection of the initial values of n and m as will be shown 
later . 

If in equation (2) the initial conditions are set equal 
to zero, after a simple manipulation, the transfer function 
of the system can be obtained and is given in equation (15)-. 



U(s) 



Y(S) 




n-1 



. + b. s + b 
1 



. + a-, s + a 
1 < 



o 



(15) 



o 
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It should be noticed that the coefficients and 
of equation (10) are exactly the same as the a^ and b^ of 
equation (5). Therefore, as soon as the set of overdetermined 
algebraic equations is solved for a^ and b ^ , the transfer 
function is determined. 

Equation (15) is usually written in the form (16) . 



Y(s) _ K[s ra + c in . 1 s m_1 + ... + OjS + cj (1 6) 

U(s) s n + a . s ra + ... + a.s + a 
n-1 1 o 

Id 

Where K = b is the so called constant gain and c = g for 

m g ET 

m 

o < g < m. 

If the roots of the numerator and denominator polynomials 
are found, the more familiar form of the transfer function 
can be written; equation (17) . 



Y(s) = 
U(s) 



K 



m 

* = l x (s + q j ) 

h (s + p. ) 

nq 



(17) 



Where -p . and 
Equation 
form as 



- q j 

(16) 



, are the poles and zeros respectively, 
can also very easily be written in state 
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y 




* (17a) 

! x 

L. nl 

But, how are n and m determined? It was said that in 
order to apply the method efficiently a rough idea of the 
order of the system should be available, although it is not 
absolutely necessary; that rough idea is used to pick n' and 
m 1 as the first approximation for n and m, taking into account 

that n'>n and m'>m to be sure that neither poles nor zeros 
are missing. Since m is usually smaller than n, m' should 
not be greater than n'-l. 

Consider then that n and m are the true numbers of poles 
and zeros and that n' and m' are assumed. Under these 
conditions one wbuld expect to obtain n'-n, or m'-m - what- 
ever number is smaller - poles equal to zeros, what is 
equivalent to saying that they are not present in the system, 
what is actually the case; the (n' - n) - (m' - m) poles or 
(m'- n) - (n' - n) zeros in excess which have not been 
cancelled will appear to be either in the right half plane 
or with values too big compared with the actual values of the 
poles and zeros in the system. As a result after a first 
guessing for n and m, a clear indication of the actual order 
of the system and of the number of zeros, should show up. 
Examples in Chapter III will demonstrate this point. 
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B. IDENTIFICATION OF THE INITIAL CONDITIONS 

Although the knowledge of the initial conditions is 
not important by itself, it might be if one desired to 
compare the system model with the actual system by exciting 
the system model with the same input that was used in the 
identification. By comparing the output of the model with 
the output of the system a rough idea of the accuracy of 
the model can be obtained. However, in order for this to be 
possible, all the initial conditions of the input and output 
should be known. 

When the overdetermined set of linear algebraic equations 
is solved, besides the values of a. and b., the d, initial 
condition terms are also found. These d^ terms are related 
with the initial conditions by equation (18) . This equation 
can easily be found by specifically writing in the contribu- 
tion of the individual initial conditions when the Leplace 
Transform of equation (1) was taken, to obtain equation (2) . 



"d 

o 




a l a 2 


. ~\ 

• • • cl ^ X 

n-l ^ 




“ Y (o) 


d l 




a 2 a 3 


/ 




y l)(o) 
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m — 
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v / 
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1 • 
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• 
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/ A 
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• 








1 J 




n-2 ) , , 
y (o) 


d 0 

n-Z 




a 1 ^ 
n-l x 


\J 




d i 

n-l 




1 ' 






y (o) 


- 




* 






» — 



( 18 ) 
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1 0 ... 0 




u (o) 


• • 
I • • 




. u 1 ^ (o) 
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• • 

1 • • 

• • 
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• 

m-1 , N 
u (o) 


. 0 ... 0 






• 0 ... 0 
1 • • 
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1 # • 

| 0 ... 0^ 
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Equation (18) requires the knowledge of m-1 derivatives 
of the input and numerical techniques may be necessary to 
find them, however, this difficulty may be avoided in many 
cases by choosing the beginning of the input-output record 
at a moment where the input is constant and thus all the 
derivatives would be zero. 



C. SPECIAL CASES WHERE SIMPLIFICATION COULD BE POSSIBLE 

Theoretically the method should identify the system 
independently of the number of unknowns, but practically 
this is not the case due to the natural limitations in the 
accuracy of the numerical integration methods and in the 
measure of the input and output records, thus it would be 
desirable to reduce the number of unknowns as much as 
possible, because as a general rule the fewer the unknowns, 
the more accurate its identification. 

A great reduction in the number of unknowns occurs if 
the record of the input and output begins when all the 
initial conditions are zero. Under these conditions, all 



( 18 ) 



22 



the terms d^ of equation (5) are zero and the number of 
unknowns is reduced' from 2n+m+l to n+m+1 since it was 
assumed n>m. 

Any additional information could be used to simplify 
the problem, for instance if the value of the steady state 
gain constant were known the number of unknowns could be 
reduced by one. 
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'III. IMPLEMENTATION 



A. NUMERICAL METHODS 

As was mentioned before, after having the record of the 
input and output of the system, which can be obtained as 
described in Appendix I, it is necessary to perform the 
multiple integrations indicated in equation (5) , for 
different times t^, in order to form the set of overdetermined 
linear algebraic equations whose solution gives the values 
of the parameters a^, b^ and d^, although only the a^ and b^ 
are of interest in this work. The d^, even available are 
not used. 

The multiple integrations can be handled for any numerical 
method subroutine and although there are quite a few highly 
sophisticated numerical integration procedures available, 
trapezoidal integration will be used because most of the 
more complex integration techniques perform poorly when the 
function being integrated is discontinuous. Since the 
inputs used here are step-wise waves, complex integration 
techniques are not desirable. On the other hand, the slight 
improvement in accuracy offered by more advanced methods is 

not, in general, enough to justify the tremendous increase 
in computational load associated with their use. 

The solution of the set of overdetermined linear equations 
is a classical problem in several fields of mathematics and 
engineering. Most of the classical techniques for solving 
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such problem are not practical. They tend to magnify 
the errors introduced by the finite precision of the computer 
to the point where the solution is meaningless. Fortunately, 
several modern methods are available that display more 
acceptable behaviour. The method used here was developed 
by Golub [22] . The basic approach is to triangularize the 
coefficient matrix by performing a Choleski decomposition. 

The decomposition is accomplished by applying repeated 
Householder transformations [23] . Once the coefficient 
matrix has been triangularized the unknowns can be obtained 
by back substitution. The method is quite stable numerically 
and is capable of handling ill-conditioned coefficient 
matrices. 

B. CHARACTERISTICS OF GOOD . INPUT-OUTPUT RECORDS 

The accuracy with which a system can be identified is 
strongly dependent on the input-output record used in the 
identification. Since parameters are identified on the 
basis of their effect on the output, it will be impossible 
to identify a parameter unless its effect is measurable. If 
the input to a system has a frequency spectrum that is more 
or less uniform over the frequency range of interest, the 
identification will probably be very good. If the frequency 
spectrum of the input is confined to a narrow band of fre- 
quencies, the identification will probably be very bad. It 
is well known that signals with sharp discontinuities have 
a broader bandwidth than slowly varying continuous signals. 
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For this reason, the step-wise wave was chosen as the input 
in this research, after having demonstrated experimentally 
the preceeding observations seeing that with a step input 
it was never possible to obtain a good identification even 
for a second order system. Periodic functions should also 
be avoided at the input, because very probably many of the 
algebraic equations are going to be linearly dependent. This 
fact was also proved experimentally; a periodic square wave 
was tried as input and the identification was also bad when 
a step-wise input gave very good results. 

Another desired characteristics of the input, experimentally 
proved, is that it should be adjusted in such a way that the 
output, if it is strictly increasing, should at some moment 
in the sampling interval have a negative second derivative; 
in other words, -the output is desired to have an inflexion 
point within the period in which the record is taken. The 
reason for this is thought to lie in the subroutine DLLSQ 
which solves the overdetermined set of algebraic equations. 

What DLLSQ essentially does is look for a fitting output 
function that minimizes the least square error between the 
generated and the measured and if the second and higher 
order derivatives of the output are positive at the end of 
the sampling interval, the function generated by the sub- 
routine may in some cases correspond to an unstable situation. 
Examples in Chapter IV will also demonstrate this point. 

In relation to the sampling rate at which the input and output 
should be digitized, it can be said that low sampling rates 
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may introduce inaccuracies in the numerical integrations 
and as a result poor identification, especially in the 
location of high frequency poles and zeros. Maximum 
sampling rate available or ten to one hundred times shorter 
than the shortest system time constant is recommended. The 
output of the system is also recommended to have maximum 
value of the order of the maximum value of the input, thus 
the numbers with which the computer deals are not very 
different in magnitude. If the output has not the ' same order 
of magnitude of the input, it can as a rule be amplified 
accordingly, taking this amplification factor into account 
when the time comes to see what the actual identification is. 

C. PROGRAMS DESCRIPTION 

Four different computer programs were used in this 
research. A listing of them is given after Chapter IV. The 
first one is the A/D program. This program is used for high 
speed A/D (analog to digital) conversion. Its output is a 
seven track magnetic tape that contains the digitized input 
and output record of the system. This program consists of 
a main driver program and five subroutines. 

The source language is FORTRAN for the main program and 
METASYMBOL for the subroutines. 

The A/D program provides a means of sampling and digitizing 
analog information at a rate (up to lOKHz) determined by an 
interrupt signal originated at the analog computer logic 
patchboard (Ci 5000) . (See Appendix I) . The digital data 
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are stored sequentially in one of two buffers provided by 
the main program. When one buffer is full, the storage 
sequence is continued at the top of the other buffer. An 
indicator notifies the main program of this switching and 
the main program processes the data in the full buffer 
writing them in the magnetic tape while the other buffer 
is being filled. 

The parameters to be entered in the main program are: 

NS AMP = Number of words intone block/NCHAN . 

NCHAN = Number of channels to be digitized. 

NREC = Maximum number of records to be digitized in 
one file. 

ITAPE = Number of tape drive unit to be used. 

NDEL = Delay between digital and analog signal in 
option 7. The options of this program are: 

1 = Desired to input new parameters. 

2 = Start digitizing the analog input signals. 

3 = Write an END OF FILE on the magnetic tape. 

4 = Rewind the tape to the loading point (beginning) . 

5 = Skip the following number of END OF FILE marks. 

6 = Print the following amount of digitized data. 

7 = Actuate the A/D program for the next single block 
of data encountered and feed the signal to the recorder. 

The five subroutines are: ADSTART, ADSTOP, ADFAST, 

MTRDY and MTOUT. 

Subroutine ADSTART initializes the A/D input operations. 
The calling sequence is: 
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CALL ADSTART (NCH, BUF , NEWBUF, NSAM, RECNUM, NEXLOC) , 
where 

NCH = Number of A/D channels to sample. 

BUF = A pointer to the first buffer. 

NEWBUF = A pointer to the second buffer. 

NSAM = Number of samples per buffer. 

RECNUM = A record counter which is bumped by ADFAST 
each time a buffer is filled. 

NEXLOC = The return address for ADFAST to be used each 
time a buffer is filled. 

The first four arguments are used to set up pointers 
and counters in ADFAST. RECNUM provides a means for ADFAST 
to keep the main program aware of the record number of the 
buffer currently being filled. NEXLOC is an alternate 
return location for ADFAST. ADSTART also sets up the 
interrupt controller branch instructions and starts the A/D 
hardware. ADSTOP stops the A/D operations and restores the 
interrupt locations to their normal instructions. 

ADFAST controls the buffer location pointer which directs 
the A/D data into core. 

MTRDY initializes the magnetic tape handler. It has 
the following arguments : 

UNIT = Magnetic tape unit number. 

BUFl = First word address of first buffer. 

BUF2 = First word address of second buffer. 

NW = Number of words to be written per record. 

IR = Tape operation status word return. 
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1 if tape I/O channel is busy. 

2 if operation completed. 

3 if tape unit is not ready. 

4 if tape writes error. 

MTOUT controls the stream of A/D data from the buffers 
into the tape. It is called by the main program each time 
a record of data is ready for recording. The argument, NB, 
indicates which buffer is just filled and is to be output. 
After setting the status indicator ^IR) to 1, to indicate 
an operation not complete condition, the routine tests the 
tape and channel status. If either is busy, IR is bumped to 
3 (not ready condition) and the routine exits back to the 
calling program. If everything is all right, NB is checked 
and the I/O instructions for writing the filled buffer on 
tape are executed. If the operation is completed without 
error, IR is set to 2 (operation completed) and to 4 (write 
error) if an error occurred, and execution is returned to 
main program. 

The seven track magnetic tape generated by the A/D 
program recording the digitized input and output can not be 
directly read by the IBM 360 computer system, and since this 
system is used to perform the integrations of equation (5) , 
a seven to nine track tape conversion should be done prior 
to the numerical integration. This conversion is done by 
the IBM 360 using the CONVERT program, which is described 
in Ref. 21. This program has several parameters. One of 
them is BLKSIZE which is related with the NSAMP and NCHAN 
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parameters used in the A/D program. BLKSIZE can not be 
greater than 32.700 and for NCHAN = 2 limits the maximum 
number of samples that can be taken for the input and out- 
put record to about 4000. 

The third program is the main one and is a direct 
implementation of the procedure developed in Chapter II. 

► 

This program has several important parameters that have to 
be defined at the very beginning. NP and NZ are rough 
estimates of the number of poles and zeros in the system to 
be identified. NP and NZ stand for n' and m' in the program, 
N' and m' were defined in Chapter II. KPTMAX is the number 
of sample points available in the input-output record. The 
input amplitude is R and the output amplitude is C. Because 
the integration of the time has to be done, the program 
generates the time corresponding to each sample for what 
the sampling rate has to be known and is represented by the 
parameter DT in the program. IPTS is the number of sample 
points that will separate successive linear equations. 
Therefore, the total number of linear algebraic equations 
formed is equal to KPTMAX/IPTS. The program reads- in the 
nine track tape generated by the CONVERT program, performs 
the integrations of equation (5) and forms the overdetermined 
set of linear algebraic equations. When the set is formed, 
subroutine DLLSQ is called. This subroutine is an implemen- 
tation of the Golub algorithm for solving overdetermined 
sets of linear equations. Subroutine DLLSQ returns the 

values a., b. and d, . The a. and b. coefficients are fed 
i n k i i 
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into the RTPLSB subroutine, which is a polynomial root 
finder, in order to obtain the poles and zeros of the system. 

The output of the program are the results of the 
identification given in both transfer function and state 
variable form. 

The fourth program was written to compare the outputs 
of the actual system and of the ones obtained from the 
identification, when all of them are excited with a step. 

This program solves the differential equations that describe 
the system. The differentials equations should be arranged 
in state variable form. The output of the program is a 
drawing of the response of the different systems to be 
compared in the same plot. 

D. EXAMPLES 

The following examples show the possibilities of the 
method, how the order of the system can in general be 
determined, if it is not already known, from a trial identi- 
fication run using an estimated order greater than the 
actual order of the system. They also show how the accuracy 
of an identification depends on the input wave form, on the 
R.M.S. value of the noise present, and as expected, does 
not depend on the initial conditions. They also indicate 
how the method can be used to reduce or obtain a lower order 
mathematical model. In all the examples, the MAIN PROGRAM 
PARAMETERS were: 

NUMBER OF DATA POINTS = 4000 
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NUMBER OF EQUATIONS = 40 

DATA POINTS PER EQUATION = 100 

In examples 1 through 6, the input function was a step- 
wise increasing wave form (staircase wave form, see Fig. 1, 
Example 1), with a maximum value of about 85 volts. The 
amplitude and the number of steps were changed in each 
example. In Example 3 and 4, the maximum amplitude of the 
input was about 9 volts. In examples 6 through 10, the 
input function was also a step-wise wave form, but instead 
of being strictly increasing, it had one part strictly 
increasing and the final part strictly decreasing. See 
Fig. 4, Example 5. 

The systems investigated in Examples 1, 2, 5, 6, 7, 8, 

9, 10 and 11 were set up using R-C and R-C-L circuits 
isolated by operational amplifiers which also provided the 
desired amplification factor. In Examples 3 and 4, an 
active system was studied. This system was a D.Ci positioning 
servo. 

Example 1 was performed as a first test and in it a 
second order system - two real poles - was considered. The 
input was a step-wise increasing function, (see Fig. 1) , 
with zero initial conditions. The results of the identifica- 
tion were within the 2% for the dominant pole and the gain 
constant and within the 10% for the higher frequency pole. 

This example shows how assuming different values for m' and 
n', the results of the identification gave a clear indication 
of what the actual order of the system was . 
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Example 2 illustrates the identification of a third 
order system. The input function as in Example 1 was a 
step-wise increasing function. The. initial conditions were 
also zero. The identification of the low frequency poles 
and the gain constant were also within 3% of the true values; 
however, the pole at -100.0, was identified with an error of 
about 30%, this result was not a surprise since the more 
negative the value of a pole, compared with the other poles 
in the system, the less exact its identification. This 
example shows also the fact that in some cases the trial 
runs to find out the true values of n' and m' , did not indicate 
very clearly what these values were. In cases like this, a 
thorough search should be done around the expected values of 
n and m and different input-output records are recommended 
to be used. 

This example also shows how the method could be used to 
reduce the order of the system, i.e., when n' = 3 and m' = 

1 were assumed, the results correspond practically to a 
second order system, since the pole at -2.22 and the zero at 
-2.20 cancel each other. When n' = 2 and m' = 0 were 
assumed, practically the same results were obtained. The 
steady state error coefficient of the second order has 
practically the same value of the third order. Essentially, 
the method neglected the effects of the higher frequency 
pole. 

The experiment described in Example 3 was performed 
using a D.C. positioning servo. The transfer function of 
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the servo was first obtained from the asymptotic approxima- 
tion of the closed loop frequency response plot shown in 
Fig. 2. The parameters of the servo were adjusted in such 
a way that two dominant complex roots were present. During 
the frequency response test, it was found that the servo 
had nonlinearities since a jump in the resonance frequency 
was observed and also that reducing the amplitude of the 
input sine wave from twenty volts to one volt peak to peak 
the resonance frequency was shifted from 2.6Hz to 8.1Hz 
(52 rad/sec) . 

In taking the input-output record in this example, the 
input was adjusted to be a step-wise increasing function, 

(see Fig. 1), but with a step amplitude of about 0.5 volts. 
The 4 volts step amplitude used in Examples 1 and 2 was too 
big; the D.C. motor ran continuously. 

The results of the identification were almost identical 
to the ones obtained from the frequency response plot. 

In Example 4, the same D.C. servo was treated but now 
the parameters of the servo were adjusted in order to have 
two real poles. Also a frequency response test was performed 
and a transfer function was obtained from the frequency 
response curve asymptotic approximation (see Fig. 3) . The 
results of the identification were very close to the ones 
found from the frequency response plot. 

In Example 5, a third order system is again considered, 
but now a zero was added. 
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The input function and the I.C. were the same as in 
Example 1. The results in this case also showed a poor 
identification of the higher frequency pole. 

Example 6 shows the effect of the input function in 
the accuracy of the identification. Now a fourth order 
system is studied using two different inputs. The first 
input was the same as in Example 1 and is shown together 
with the output of the system in Fig. 4. With this input, 
the results of the identification were very poor since not 
a single pole was correctly identified. Several fourth 
order systems were treated with this type of input and the 
identification was bad, giving in some cases poles in the 
right half plane. It was observed, however, that the output 
of the system, at the end of the record had a positive 
second derivative (see Fig. 4) . It was thought that perhaps 
this fact could be the cause of the poor identification 
considering that an unstable system might have the same 
output, for the finite interval of time considered, and this 
unstable system could be the one identified or generated 
by the DLLSQ subroutine. With the preceding in mind, the 
input function was changed and adjusted in such a way that 
the output had the second derivative negative at a certain 
instant within the time interval considered (see Fig. 5). 

The results in this case were very good since even the high 
frequency pole was identified with an error within the 6% 
of the true value. The I.C. were also zero value. 
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In Example 7, the identification of a system with four 
poles and two zeroes is presented. The input function was 
as the one shown in Fig. 5. The results in this case were 
also very good. 

Example 8 shows how the I.C. does not affect the identi- 
fication, the results being practically the same with I.C. 
equal to zero as with I.C. 0; several different initial 
conditions were considered. 

In Example 9, a third order system was studied, but noise 
was added at the input and output. The results showed the 
great influence of the noise in the identification. It was 
found that in order to have a good identification, the 
maximum R.M.S. value of the nois'Gf should not be greater than 
about .5 volts. 1.0 volts R.M.S. noise made the identifica- 
tion very bad. The noise added was a gaussian, low frequency 
noise - from a few cycles per second to about 300 cycles per 
second - withuzero mean. 

Examples 10 and 11 show how the method can be used to 
reduce the order of the system. In Example 10, a fifth 
order system was treated. The identification was not good 
for the higher frequency poles. Other fifth order systems 
were studied and in some cases none of the poles was identi- 
fied, which indicates that fourth order systems are perhaps 
the highest order systems that can be identified with the 
kind of inputs used in this work. 

In Example 10, the reduction was not, as a direct conse- 
quence of the poor identification, very good as can be seen 
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comparing the responses to a step input of the different 
models, see Fig. 6. However, it is observed that the 
three models generated by the method, 5^, 3 r< ^, and 2 n< ^ 
order, give responses fairly close to one another. 

In Example 11, although the identification is considered 
very good, the responses of the different reduced models 
practically coincide during the transient period, but in 
steady state an error of about 30% is observed between the 
step responses of the fourth order identified model and the 
third and second order reduced models (see Fig. 7). However, 
if the gain constant of the reduced models is adjusted in 
order to have the same steady state gain as the fourth order 
model, the responses are practically the same (see Fig. 8). 

It is to be noticed that the systems considered true, 
in the examples, were true within the tolerances of the values 
of the resistors, capacitors, coils and amplifiers used. 
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EXAMPLE I 



SYSTEM TO BE IDENTIFIED 
SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 


1 


o 

• 

o 

cH 

1 


0.0 


2 


-100.0 


0.0 


ZEROES 


REAL 


IMAGINARY 


GAIN CONSTANT = 


1000.0 





TESTS FOR EXAMPLE I 

NO. OF ZEROES RESULTS 

0 Good 

1 See Table of Data 

2 See Table of Data 

1 See Table of Data 

INPUT FUNCTION STEP-WISE STRICTLY INCREASING SHOWN IN 
FIG. 1. 

INITIAL CONDITIONS EQUAL ZERO. 



RUN 

a 

b 

c 

d 



NO. OF POLES 
2 

3 

4 
2 
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INPUT 

VOLTS 




I2TPUT FUITCTIOIT USED III EXAMPLE ITO^ 1 



OUTPUT 

VOLTS 




OUTPUT FUNCTION EXAMPLE 110. 1. 
FIGURE 1 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 



SYSTEM TRANSFER FUNCTION 
POLES REAL 

• 1 -9.8488868 

2 -111.04012 



IMAGINARY 

0.0 J 

0.0 J 



ZEROES 



REAL IMAGINARY 



GAIN CONSTANT = 984.74146 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 

1 1093.6215 

2 120.88900 

B VECTOR 

2 984.74146 



C VECTOR 

1 1.0000000 



REMARKS : 

ASK FOR TWO POLES NO ZEROES. 
GOOD RESULTS. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 
SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




. 1 


-97.994740 


o 

• 

o 


J 


2 


2.2362954 


o 

• 

o 


J 


3 


-9.8646060 


o 

• 

o 


J 


ZEROES 


REAL 


IMAGINARY 




1 


2.2427425 


o 

• 

o 


J 



GAIN CONSTANT =• 867.90796 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 

1 -2161.7809 

2 725.47414 

3 105.62305 

B VECTOR 

3 867.90796 



C VECTOR 

1 -2.2427425 

2 1.0000000 



REMARKS : 

ASK FOR THREE POLES ONE ZERO. 

RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE SYSTEM 
SINCE ZERO AT 2.24 CANCELS POLE AT 2.23. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (c) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


-4.4110354 


-54.484437 


J 


2 


-4.4110354 


54.484437 


J 


3 


-99.827986 


o 

• 

o 


J 


4 


-9.8612480 


o 

# 

o 


J 


ZEROES 


REAL 


IMAGINARY 




1 


-7.0865674 


-59.922949 


J 


2 


-7.0865674 


59.922949 


J 


GAIN CONSTANT 


= 727.27368 






SYSTEM STATE VARIABLES (PHASE 


FORM) 




A VECTOR 








1 


2941483.4 






2 


336437.35 






3 


4940.1259 






4 


118.51130 






B VECTOR 








4 


727.27368 






C VECTOR 








1 


3640.9793 






2 


14.173135 






3 


1.0000000 
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REMARKS : 



ASK FOR 4 POLES TWO ZEROES. 

RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE 
SYSTEM, SINCE COMPLEZ ZEROES CANCEL COMPLEX POLES. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (d) 



SYSTEM TRANSFER FUNCTION 
POLES REAL 

1 -9.8592476 

2 -97.619069 



IMAGINARY 

0.0 

0.0 



ZEROES 

1 



REAL 

-847.73348 



IMAGINARY 

0.0 



GAIN CONSTANT = 1.0221453 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 



1 

2 

B VECTOR 
2 

C VECTOR 
1 
2 



962.47832 

107.47832 

1.0221463 

847.73348 

1.0000000 



J 

J 



J 



REMARKS : 

ASK FOR TWO POLES ONE ZERO. 

RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE SYSTEM. 
SINCE THE ZERO AT - 847.7 CAN BE NEGLECTED. 
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EXAMPLE 2 



SYSTEM TO BE IDENTIFIED 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 




IMAGINARY 






I 


-1.000 




o 

• 

o 




J 


2 


-10.00 




o 

• 

o 




J 


3 


-100.0 




o 

• 

o 




J 


ZEROES 












GAIN CONSTANT - 10000.0 












TESTS FOR EXAMPLE 2 








RUN NO. OF POLES 


NO. OF ZEROES RESULTS 




a 


3 


0 


Good 




b 


4 


1 


See Table 


of 


Data 


c 


4 


0 


See Table 


of 


Data 


d 


4 


2 


See Table 


of 


Data 


e 


3 


1 


See Table 


of 


Data 


f 


3 


1 


See Table 


of 


Data 



INPUT FUNCTION STEP-WISE STRICTLY INCREASED SIMILAR 
TO THE ONE SHOWN IN FIG. 1. 

INITIAL CONDITIONS EQUAL ZERO. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 
SYSTEM TRANSFER FUNCTION 



POLES 


REAL IMAGINARY 


1 


-1.1072321 0.0 J 


2 


-10.224865 0.0 J 


3 


-131.11147 0.0 J 


ZEROES 


REAL IMAGINARY 


GAIN CONSTANT 


= 11153.488 



SYSTEM STATE VARIABLES (PHASE FORM) 



A VECTOR 




1 


1484.3522 


2 


1497.0893 


3 

B VECTOR 


142.44357 


3 


11153.488 


C VECTOR 


1.0000000 



REMARKS : 

ASK FOR THREE 


POLES NO ZEROES. 



GOOD RESULTS FOR LOW FREQUENCY POLES, NOT SO GOOD FOR 
HIGH FREQUENCY POLE. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 
SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




■ 1 


-1.1675351 


o 

• 

o 


J 


2 


-9.9863699 


o 

• 

o 


J 


3 


-1.7691065 


-59.923677 


J 


ZEROES 


REAL 


IMAGINARY 




1 


157.78387 


o 

• 

o 


J 



GAIN CONSTANT = -1919.7240 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 



1 

2 

3 

4 

B VECTOR 
4 

C VECTOR 
1 
2 



41903.747 

40128.130 

3645.1011 

14.692118 

-1919.1240 

-157.78387 

1.0000000 



REMARKS : 

ASK FOR FOUR POLES ONE ZERO. 

RESULTS TO NOT GIVE CLEAR INDICATION OF THE ORDER OF THE 
SYSTEM, INDICATING FOUR POLES WHEN THERE ARE THREE. THE 
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ZERO IN THE RIGHT HALF PLANE SHOULD BE NEGLECTED. 



49 










4 



* 
















IDENTIFICATION OF UNKNOWN SYSTEM - Run (c) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


-1.1300231 


o 

• 

o 


J 


2 


-10.158884 


0.0 


J 


3 


-11.337775 


-53.285012 


J 


4 


-11.337775 


+53.285012 


J 


ZEROES 


REAL 


IMAGINARY 





GAIN CONSTANT = 252014.00 



SYSTEM STATE VARIABLE (PHASE FORM) 
A VECTOR 



1 

2 

3 

4 

B VECTOR 
4 

C VECTOR 

1 



34070.104 

33763.954 

3235.2996 

33.964457 

252014.00 

1.0000000 



REMARKS : 

ASK FOR FOUR POLES NO ZEROES. 

RESULTS DO NOT GIVE CLEAR INDICATION OF ORDER OF THE 
SYSTEM INDICATING FOUR POLES WHEN THERE ARE THREE. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (d) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


-1.1196795 


o 

• 

o 


J 


2 


-10.190750 


o 

• 

o 


J 


3 


0.16021282 


-59.917890 


J 


4 


0.16021282 


59.917890 


J 


ZEROES 


REAL 


IMAGINARY 




1 


3.5557973 


-30.108372 


J 


2 


3.5557973 


30.108372 


J 


GAIN CONSTANT = 332.37549 






SYSTEM STATE 


VARIABLES (PHASE 


FORM) 




A VECTOR 








1 


40965.288 






2 


40602.813 






3 


3597.9654 






4 


10.990004 






B VECTOR 








4 


332.37549 






C VECTOR 








1 


919.15775 






2 


-7.1115946 






3 


1.00000000 
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REMARKS : 



ASK FOR FOUR POLES TWO ZEROES 

RESULTS INDICATE TWO POLES WHEN THERE ARE THREE. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (e) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


-1.0459484 


o 

• 

o 


J 


2 


-10.251872 


o 

• 

o 


J 


3 


2.2244020 


o 

• 

o 


J 


ZEROES 


REAL 


IMAGINARY 




1 


2.2033545 


o 

• 

o 


J 



GAIN CONSTANT = 83.810272 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 



1 

2 

3 

B VECTOR 
3 



23.852104 

35.853822 

13.522222 

83.810272 



C VECTOR 



1 2.2033545 

2 1.0000000 



REMARKS : 

ASK FOR THREE POLES ONE ZERO. 

RESULTS INDICATE TWO POLES WHEN THERE ARE THREE. 
POLE AT -2.2 IS CANCELLED BY ZERO AT -2.22. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (f) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 


1 


-1.0934113 


o 

• 

o 


2 


-10.207312 


o 

• 

o 


ZEROES 


REAL 


IMAGINARY 


GAIN CONSTANT 


= 90.3204350 


SYSTEM STATE VARIABLES (PHASE 


FORM) 


A VECTOR 


1 


11.232745 




2 


11.300723 




B VECTOR 


• 




2 


90.320435 




C VECTOR 


1 


1.0000000 





REMARKS : 

ASK FOR TWO POLES NO ZEROES. 

RESULTS INDICATE TWO POLES NO ZEROES WHEN THERE ARE 
THREE. METHOD NEGLECTED POLE AT -100.0, OTHERWISE IDENTI- 
FICATION VERY GOOD. 
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EXAMPLE 3 



SYSTEM TO BE IDENTIFIED 
SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


1 

• 

00 


-51.5 


J 


2 


1 

• 

00 


51.5 


J 


ZEROES 


REAL 


IMAGINARY 




GAIN CONSTANT = 
RUN NO . OF 


= 2950.0 

TESTS : 
POLES 


FOR EXAMPLE 3 

NO. OF RESULTS 


RESULTS 


a 2 

D.C. POSITIONING SERVO 


0 

FREQUENCY RESPONSE 


Good 

IN FIG. 2 



INPUT FUNCTION STEP-WISE INCREASING SIMILAR TO THE ONE 
SHOWN IN FIG. 1, BUT WITH STEP AMPLITUDE ABOUT Q.5 VOLTS. 
INITIAL CONDITIONS EQUAL ZERO. 
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Frequency Response System Example 
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10.0 52.0 100.0 ioio .0 



IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL IMAGINARY 


1 


-5.3585933 -56.771298 J 


2 


-5.3585933 56.771298 J 


ZEROES 


REAL IMAGINARY 


GAIN CONSTANT 


= 3332.2095 



SYSTEM STATE VARIABLES (PHASE FORM) 



A VECTOR 




1 


3251.6948 


2 

B VECTOR 


10.717187 


2 

C VECTOR 


3332.2095 


1 


1.0000000 


REMARKS: 





ASK FOR TWO POLES NO ZEROES. 
GOOD RESULTS. 
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EXAMPLE 4 



SYSTEM TO BE IDENTIFIED 
SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


1 

• 

VO 


o 

• 

o 


J 


2 


-32.0 


o 

• 

o 


J 


ZEROES 


REAL 


IMAGINARY 




GAIN CONSTANT = 


150.0 







TESTS FOR EXAMPLE 4 

RUN NO. OF POLES NO. OF ZEROES RESULTS 

a 2 0 Good 

D.C. POSITIONED SERVO. 

FREQUENCY RESPONSE IN FIG. 3. 

INPUT FUNCTION THE SAME AS IN EXAMPLE 2. 

INITIAL CONDITIONS EQUAL ZERO. 



58 



Frequency Response System Example 
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100 



IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 



SYSTEM TRANSFER FUNCTION 
POLES REAL 

. 1 -5.2129680 

2 -36.525846 



IMAGINARY 

0.0 

0.0 



J 

J 



ZEROES 



REAL 



IMAGINARY 



GAIN CONSTANT = 178.37120 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 



1 

2 

B VECTOR 

2 

C VECTOR 

1 



190.40807 

41.738815 

178.37120 

1.0000000 



REMARKS : 

ASK FOR TWO POLES NO ZEROES. 
GOOD RESULTS. 
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; 




























EXAMPLE 5 





SYSTEM TO : 


BE IDENTIFIED 




SYSTEM TRANSFER 


FUNCTION 






POLES 


REAL 


IMAGINARY 




1 


-1.0 


0.0 


J 


2 


-10.0 


0.0 


J 


3 


-55.0 


0.0 


J 


ZEROES 


REAL 


IMAGINARY 




1 


' -5.0 


0.0 


J 


GAIN CONSTANT = 500.0 








TESTS FOR 


EXAMPLE 5 




RUN NO . 


OF POLES 


NO. OF ZEROES 


RESULTS 



a 3 1 Good 

b 2 0 See Table of Data 

INPUT FUNCTION STEP-WISE STRICTLY INCREASING SIMILAR TO 
THE ONE SHOWN IN FIG. 1. 

INITIAL CONDITIONS EQUAL ZERO. 
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< 







IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 


1 


-0.89534766 


o 

• 

o 


2 


-8.3914174 


o 

• 

o 


3 


-88.600646 


o 

• 

o 


ZEROES 


REAL 


IMAGINARY 


1 


-3.8670622 


o 

• 

o 


GAIN CONSTANT 


= 706.88770 




SYSTEM STATE VARIABLES (PHASE 


FORM) 


A VECTOR 


1 


665.67756' 




2 


830.32662 




3 


97.887411 




B VECTOR 


3 


706.88770 




C VECTOR 


1 


3.8670622 




2 


1.0000000 





J 

J 

J 



J 



REMARKS : 

ASK FOR THREE POLES ONE ZERO. 

RESULTS SHOW POOR IDENTIFICATION OF THE HIGHER 
FREQUENCY POLE. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 



SYSTEM TRANSFER FUNCTION 
POLES REAL 

1 -3.2441577 

2 -108.41738 



IMAGINARY 

0.0 

0.0 



ZEROES 



REAL 



IMAGINARY 



GAIN CONSTANT = 763.43213 



J 

J 



SYSTEM STATE VARIABLES (PHASE FORM) 

A VECTOR 

1 351.72309 

2 111.66154 
B VECTOR 

2 763.43213 

C VECTOR 

1 1.0000000 

REMARKS : 

ASK FOR TWO POLES NO ZEROES. 

RESULTS CAN BE CONSIDERED AS A LOWER ORDER MATHEMATICAL 
MODEL. 



63 



I 







EXAMPLE 6 



SYSTEM TO BE IDENTIFIED 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


o 

• 

CM 

1 


o 

• 

o 


J 


2 


O 
• 
o 
1 — 1 

1 


o 

• 

o 


J 


3 


-33.0 


o 

• 

o 


J 


4 


-100.0 


o 

• 

o 


J 


ZEROES 


REAL 


IMAGINARY 




GAIN CONSTANT = 


= 660000 








TESTS 


FOR EXAMPLE 6 




RUN NO. OF 


POLES 


NO. OF ZEROES 


RESULTS 



a 

b 

c 

d 

e 

f 

g 

h 



4 

4 

6 

4 

5 

4 

5 
3 



0 Bad 

0 Good 

2 See Table of 
2 See Table of 
2 See Table of 

1 See Table of 
1 See Table of 
0 See Table of 



Data 

Data 

Data 

Data 

Data 

Data 



INPUT FUNCTION FOR RUN (a) STEP-WISE STRICTLY INCREASING 
SHOWN IN FIG. 4. 

INPUT FUNCTION FOR RUNS (b) , (c) , (d) , (e), (f), (g) , and 
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(h) STEP-WISE INCREASING-DECREASING SHOWN IN FIG. 5. 
INITIAL CONCLUSION EQUAL ZERO. 
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INPUT 

VOLTS 




FIRST INPUT FUNCTION USED IN EXAMPLE NO. 6 



OUTPUT 

VOLTS 




OUTPUT C ORRE ~ PON-lENT TO THE ANO /E INPUT- 
FIGURE 4 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


-1.5013900 


-4.3496191 


J 


2 


-1.5013900 


4.3496191 


J 


3 


-5.4821332 


-26.089080 


J 


4 


-5.4321332 


26.089080 


J 


ZEROES 


REAL 


IMAGINARY 





GAIN CONSTANT = 62546.430 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 



1 

2 

3 

4 

B VECTOR 
4 

C VECTOR 
1 



15047.775 

2366.2077 

764.79050 

13.967046 

62546.430 

1.0000000 



REMARKS: 

INPUT FUNCTION - STEP-WISE INCREASING FUNCTION SHOWS IN 
FIG. 4. 

ASK FOR FOUR POLES NO ZEROES. 
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RESULTS NOT GOOD. 



NOTE THE POSITIVE SECOND DERIVATIVE OF THE OUTPUT AT 
THE END OF THE SAMPLING INTERVAL. 



\ 
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INPUT 

VOLTS 




SECOiJD INPUT FUNCTION USED IN EXAMPLE HO. 6 



OUTPUT 

VOLTS 




OUTPUT CORNS SP CEDENT TO THE ABOVE INPUT 



FIGURE 5 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


-2.0172066 


o 

• 

o 


J 


2 


-10.764637 


o 

• 

o 


J 


3 


-33.596973 


o 

• 

o 


J 


4 


-105.65671 


o 

• 

o 


J 


ZEROES 


REAL 


IMAGINARY 





GAIN CONSTANT = 612357.94 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 



1 

2 

3 

4 

B VECTOR 
4 

C VECTOR 

1 



77080.941* 

48396.117 

5351.3789 

152.03553 

612357.94 

1.0000000 



REMARKS: 

INPUT FUNCTION STEP-WISE INCREASING-DECREASING SHOWN 
IN FIG. 5. 

ASK FOR FOUR POLES NO ZEROES. 
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RESULTS VERY GOOD. 

NOTE THE NEGATIVE SECOND DERIVATIVE DURING SOME TIME 
WITHIN THE SAMPLING INTERVAL. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (c) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 




IMAGINARY 




1 


-2.0375228 




0.0 


J 


2 


-10.649015 




0.0 


J 


3 


-34.253558 




0.0 


J 


4 


-101.99791 




0.0 


J 


5 


6.7412854 




-96.372163 


J 


6 


6.7412854 




96.372163 


J 


ZEROES 


REAL 




IMAGINARY 




1 


6.2635645 




-96.069197 


J 


2 


6.2635645 




-96.069197 ' 


J 


GAIN CONSTANT = 603327.00 








SYSTEM STATE 


VARIABLES (PHASE 


FORM) 




A VECTOR 










1 


0.7075898D 


09 






2 


0. 44024818D 


09 






3 


48381252.0 








4 


1366621.4 








5 


12569.020 








6 


135.45544 








B VECTOR 










6 


603327.00 
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C VECTOR 



1 


9268.5229 


2 


-12.527129 


3 


1.0000000 


REMARKS: 





ASK SIX POLES TWO ZEROES. 

RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE 
SYSTEM, SINCE COMPLEX POLES ARE CANCELLED BY COMPLEX ERRORS 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (d) 



SYSTEM TRANSFER FUNCTION 



POLES 

1 

2 

3 

4 

ZEROES 

1 

2 



REAL 

-2.0305074 

-10.677930 

-34.308456 

-96.072665 

REAL 

-516.60587 

751.83224 



IMAGINARY 

0.0 

0.0 

0.0 

0.0 

IMAGINARY 

0.0 

0.0 



J 

J 

J 

J 



J 

J 



GAIN CONSTANT = -1.4570827 

SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 



1 

2 

3 

4 

B VECTOR 
4 

C VECTOR 
1 
2 
3 



71464.878 

44715.215 

4974.7267 

143.08956 

-1.4570827 

-388400.95 

-235.22636 

1.0000000 
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REMARKS 



ASK FOR FOUR POLES TWO ZEROES. 

RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE 
SYSTEM, SINCE THE TWO EXCESS ZEROES APPEAR ONE IN THE 
R.H.P. AND THE OTHER IN THE L.H.P. BUT CAN BE NEGLECTED. 
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V «. 







IDENTIFICATION OF UNKNOWN SYSTEM - Run (e) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


-2.7344050 


o 

• 

o 


J 


2 


-8.8026617 


o 

• 

o 


J 


3 


-43.843968 


o 

• 

o 


J 


4 


-68.587703 


o 

• 

o 


J 


5 


8.0082667 


o 

• 

o 


J 


ZEROES 


REAL 


IMAGINARY 




1 


7.9267041 


o 

• 

o 


J 


2 


-622.31546 


o 

• 

o 


J 



GAIN CONSTANT = 777.31934 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 

1 -579657.55 

2 -227126.83 

3 2737.3550 

4 3335.5841 

5 115.96047 



B VECTOR 
5 



777.31934 
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REMARKS : 



ASK FOR FIVE POLES TWO ZEROES 

RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE 
SYSTEM, SINCE THE EXCESS POLE AT 8.0 IS CANCELLED BY THE 
ZERO AT +7.9 AND THE OTHER EXCESS ZERO AT - 622.3 CAN BE 
NEGLECTED . 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (f) 



SYSTEM TRANSFER FUNCTIONS 



POLES 


REAL 


IMAGINARY 




1 


-2.0267057 


o 

• 

o 


J 


2 


-10.704298 


o 

• 

o 


J 


3 


-34.037221 


o 

• 

o 


J 


4 


-101.07630 


o 

• 

o 


J 


ZEROES 


REAL 


IMAGINARY 




1 


-4045.2488 


o 

• 

o 


J 



GAIN CONSTANT = 146.26987 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 



1 

2 

3 

4 

B VECTOR 
4 

C VECTOR 
1 
2 



74636.682 

46730.407 

5182.1817 

147.84453 

146.26987 

4045.2488 

1.0000000 



REMARKS : 

ASK FOR FOUR POLES ONE ZERO 



RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE 

78 



SYSTEM, 



SINCE THE EXCESS ZERO ZT - 4045.2 CAN BE NEGLECTED. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (g) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


-2.1810404 


o 

• 

o 


J 


2 


-10.145443 


o 

• 

o 


J 


3 


-35.910038 


o 

• 

o 


J 


4 


-96.915087 


o 

• 

o 


J 


5 


16.571798 


o 

• 

o 


J 


ZEROES 


REAL 


IMAGINARY 




1 


16.789359 


o 

• 

o 


J 


GAIN CONSTANT 


= 580130.69 






SYSTEM STATE VARIABLES (PHASE 


FORM) 




A VECTOR 


- 






1 


-1276179.1 






2 


-682609.54 






3 


-39334.691 






4 


2734.1957 






5 


128.57981 






B VECTOR 








5 


580130.69 






C VECTOR 








1 


-16.789359 






2 


1.0000000 
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REMARKS : 



ASK FOR FIVE POLES ONE ZERO. 

RESULTS GIVE CLEAR INDICATION OF THE ORDER OF 1HE 
SYSTEM, SINCE EXCESS POLE AT 16.5 IS CANCELLED BY EXCESS 



ZERO AT 16.7. 



IDENTIFICATION OF UNKNOWN SYSTEM - Run (h) 



SYSTEM TRANSFER FUNCTION 
POLES REAL 

1 -1.0271020 

2 -15.051165 

3 -15.051165 



IMAGINARY 

0.0 

-9.0280041 

9.0280041 



J 

J 

J 



ZEROES 



REAL IMAGINARY 



GAIN CONSTANT = 4051.1528 



SYSTEM STATE VARIABLES PHASE 



A VECTOR 
1 
2 
3 

B VECTOR 
3 

C VECTOR 
1 



316.39101 

338.96060 

31.129433 

4051.1528 

1.0000000 



FORM) 



REMARKS : 

ASK FOR THREE POLES NO ZEROES. 

RESULTS CAN BE CONSIDERED AS A LOWER ORDER MATHEMATICAL 
MODEL. 
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EXAMPLE 7 





SYSTEM TO BE 


IDENTIFIED 




SYSTEM TRANSFER 


FUNCTION 






POLES 


REAL 


IMAGINARY 




1 


-25.0 


o 

• 

o 


J 


2 


I 

CT> 

• 

O 


o 

• 

o 


J 


3 


-150 


-185.3 


J 


4 


-150 


185.3 


J 


ZEROES 


REAL 


IMAGINARY 




1 


-3.1 


o 

• 

o 


J 


2 


1 

O 

• 

o 


o 

• 

o 


J 



TESTS FOR EXAMPLE 7 



RUN NO. OF POLES 

a 6 

b 4 

c 4 



NO. OF ZEROES RESULTS 

3 See Table of Data 

2 Good 

3 See Table of Data 



INPUT FUNCTION SIMILAR TO THE ONE SHOWN IN FIG. 5. 
INITIAL CONDITIONS EQUAL TO ZERO. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 




IMAGINARY 




1 


-26.773050 




0.0 


J 


2 


-72.749634 




0.0 


J 


3 


-101.14084 




173.81324 


J 


4 


-101.14084 




-173.81324 


J 


5 


960.71376 




0.0 


J 


6 


52.080470 




0.0 


J 


ZEROES 


REAL 




IMAGINARY 




1 


50.693422 




0.0 


J 


2 


-3.0909162 




0.0 


J 


3 


-9.5071748 




o 

• 

o 


J 


GAIN CONSTANT = -0.33196749E 


09 




SYSTEM STATE 


VARIABLES (PHASE 


FORM) 




A VECTOR 










1 


0.39410705D 


13 






2 


0.14131408D 


12 






3 


-0 . 12683604D 


10 






4 


-43800405. 








5 


-193111.44 








6 


-710.98986 









B VECTOR 

6 -0 . 3319 6749E 09 
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C VECTOR 



1 



-1489.6708 



2 



609.25446 



3 



-38.095331 



4 



1.0000000 



REMARKS : 

ASK FOR SIX POLES THREE ZEROES. 

RESULTS INDICATE FOUR POLES TWO ZEROES, SINCE POLE 
AT 52.0 IS CANCELLED BY ZERO AT 50.7 AND POLE AT 960.7 
CAN BE NEGLECTED. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


-26.992965 


o 

• 

o 


J 


2 


-77.607958 


o 

• 

o 


J 


3 


-126.52085 


-182.98579 


J 


4 


-126.52085 


182.98579 


J 


ZEROES 


REAL 


IMAGINARY 




1 


-3.1577742 






2 


-9.376548. 


0.0 


J 


GAIN CONSTANT = 444891.63 






SYSTEM STATE 


VARIABLES (PHASE 


FORM) 




A VECTOR 








1 


0 . 10367784D 


09 




2 


5706927.6 






3 


78054.592 






4 


357.64263 






B VECTOR 








4 


444891.63 






C VECTOR 








1 


29.609022 






2 


12.534322 






3 


1.0000000 
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REI4ARKS 



ASK FOR FOUR POLES TWO ZEROES. 
GOOD RESULTS. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (c) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




• 1 


-27.385227 


o 

• 

o 


J 


2 


-75.010363 


o 

• 

o 


J 


3 


-119 . 37880 


-181.29378 


J 


4 


-119.37880 


181.29378 


J 


ZEROES 


REAL 


IMAGINARY 




1 


-3.1367869 


o 

• 

o 


J 


2 


-9.4693811 


o 

• 

o 


J 


3 


-5036.4256 


o 

• 

o 


J 


GAIN CONSTANT 


= 82.045258 






SYSTEM STATE VARIABLES (PHASE 


FORM) 




A VECTOR 








1 


96790162.0 






2 


5315200.5 






3 


73620.633 






4 


341.15318 






B VECTOR 








4 


82.045258 






C VECTOR 








1 


149599.12 






2 


63519.731 






3 


5049.0317 
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4 



1.0000000 



REMARKS : 

ASK FOR FOUR POLES THREE ZEROES. 

RESULTS GIVE CLEAR INDICATION OF THE ORDER OF THE 
SYSTEM, SINCE ZERO AT -5036.4 CAN BE NEGLECTED. 
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EXAMPLE 8 



SYSTEM TO BE IDENTIFIED 
SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


-1.0 


o 

• 

o 


J 


2 


-10.0 


o 

• 

o 


J 


3 


25.0 


o 

• 

o 


J 


ZEROES 


REAL 


IMAGINARY 





GAIN CONSTANT = 2500 



TESTS FOR EXAMPLE 8 



RUN 


NO. OF POLES 


NO. OF ZEROES 


RESULTS 


a 


3 


0 


Good 


b 


3 


0 


Good 


c 


3 


0 


Good 



INPUT FUNCTION SIMILAR TO THE ONE SHOWN IN FIG. 5. 
INITIAL CONDITIONS DIFFERENT FROM ZERO FOR THE THREE 
DIFFERENT RUNS. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 



SYSTEM TRANSFER FUNCTION 
POLES REAL 

• 1 -0.97837856 

2 -10.976897 

3 -24.694433 



IMAGINARY 

0.0 

0.0 

0.0 



ZEROES REAL IMAGINARY 

GAIN CONSTANT = 2378.9880 



SYSTEM STATE VARIABLES (PHASE 
A VECTOR 



1 

2 

3 

B VECTOR 
3 

C VECTOR 
1 



265.20737 

305.93833 

36.649709 

2378.9880 

1.0000000 



FORM) 



J 

J 

J 



REMARKS : 

ASK FOR THREE POLES NO ZEROES. 

GOOD RESULTS WITH I.C. DIFFERENT FROM ZERO. 
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* 



4 










IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL IMAGINARY 


1 


-0.96861216 0.0 J 


2 


-10.956655 0.0 J 


3 


-24.782535 0.0 J 


ZEROES 


REAL IMAGINARY 


GAIN CONSTANT 


= 2379.6726 



SYSTEM STATE VARIABLES (PHASE FORM) 



A VECTOR 




1 


263.01082 


2 


306.15108' 


3 

B VECTOR 


36.707801 


3 

C VECTOR 


2379.6726 


1 


1.0000000 


REMARKS: 

ASK FOR THREE 


POLES NO ZEROES. 



GOOD RESULTS WITH I.C. DIFFERENT FROM ZERO. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (c) 
SYSTEM TRANSFER FUNCTION 



POLES 


REAL IMAGINARY 


1 


-24.839515 0.0 J 


2 


-0.95989923 0.0 J 


3 


-10.973508 0.0 J 


ZEROES 


REAL IMAGINARY 


GAIN CONSTANT 


= 2386.2974 



SYSTEM STATE VARIABLES (PHASE FORM) 



A VECTOR 




1 


261.64609 


2 


306.95352 


3 

B VECTOR 


36.772923 


3 

C VECTOR 


2386.2974 


1 


1.0000000 


REMARKS : 

ASK FOR THREE 


POLES NO ZEROES. 



GOOD RESULTS WITH I.C. DIFFERENT FROM ZERO. 
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EXAMPLE 9 



SYSTEM TO BE IDENTIFIED 



SAME SYSTEM EXAMPLE 8. 



TESTS FOR EXAMPLE 9 



RUN 

a 

b 

c 

d 



NO. OF POLES 
3 
3 
3 
3 



NO. OF ZEROES RESULTS 



0 

0 

0 

0 



Bad 

Almost Good 
Good 
Good 



INPUT FUNCTION STEP-WISE INCREASING-DECREASING SIMILAR 
TO THE ONE SHOWN IN FIG. 5. 

INITIAL CONDITIONS EQUAL ZERO. 

LOW FREQUENCY NOISE - FROM A FEW CYCLES TO ABOUT 300. 
CYCLES PER SECOND - WAS ADDED AT THE INPUT AND OUTPUT. 
DIFFERENT R.M.S. VALUES OF NOISE WERE CONSIDERED. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 



SYSTEM TRANSFER 
POLES 
1 
2 
3 



FUNCTION 

REAL 

-0.25785674D-01 

-16.598045 

-16.598045 



IMAGINARY 

0.0 J 

-7.2418627 J 

7.2418627 J 



ZEROES REAL IMAGINARY 

GAIN CONSTANT = 2294.0054 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 

1 8.4561459 

2 328.79567 

3 33.221877 



B VECTOR 

3 2294.0054 

C VECTOR 

1 1.0000000 



REMARKS : 

ASK FOR THREE POLES NO ZEROES. 

R.M.S. VALUE OF THE NOISE 1.0 VOLTS. 
RESULTS NOT GOOD. 
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IDENTIFICATION OF UNKNOWN SYSTEM - RUN (b) 



SYSTEM TRANSFER FUNCTION 

J 
J 
J 

ZEROES REAL IMAGINARY 

GAIN CONSTANT = 2364.4692 



POLES 

1 

2 

3 



REAL 

-23.206882 

-12.305073 

-0.72464570 



IMAGINARY 

0.0 

0.0 

0.0 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 

1 206.93154 

2 311.29595' 

3 36.236600 
B VECTOR 

3 2364.4692 



C VECTOR 

1 1.0000000 



REMARKS : 

ASK FOR THREE POLES NO ZEROES. 

R.M.S. VALUE OF THE NOISE .75 VOLTS. 
RESULTS ALMOST GOOD. 



96 



IDENTIFICATION OF UNKNOV7N SYSTEM - Run (c) 



SYSTEM TRANSFER FUNCTION 

J 
J 
J 

ZEROES REAL IMAGINARY 

GAIN CONSTANT = 2289.9822 



POLES 

1 

2 

3 



REAL 

■1.1475703 

-10.808421 

■23.475558 



IMAGINARY 

0.0 

0.0 

0.0 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 



1 

2 

3 

B VECTOR 
3 



291.17727 

293.07699 

35.431549 

2289.9822 



C VECTOR 



1 1.0000000 



REMARKS : 

ASK FOR THREE POLES NO ZEROES. 

R.M.S. VALUE OF THE NOISE 0.5 VOLTS. 
RESULTS GOOD. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (d) 



SYSTEM TRANSFER FUNCTION 

J 
J 
J 

ZEROES REAL IMAGINARY 

GAIN CONSTANT = 2381.8215 

SYSTEM STATE VARIABLES (PHASE FORM) 

A VECTOR 

1 280.49226 

2 304.73267 

3 36.753190 
B VECTOR 

3 2381.8215 

C VECTOR 

1 1.0000000 

REMARKS : 

ASK FOR THREE POLES NO ZEROES. 

R.M.S. VALUE OF THE NOISE .25 VOLTS. 

RESULTS VERY GOOD. 



POLES 

1 

2 

3 



REAL 

■1.0495050 

■10.680594 

■25.023091 



IMAGINARY 

0.0 

0.0 

0.0 
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EXAMPLE 10 





SYSTEM TO 


BE IDENTIFIED 






SYSTEM TRANSFER 


FUNCTION 








POLES 


REAL 


IMAGINARY 






1 


-1.0 


0.0 




J 


2 


-10.0 


0.0 




J 


3 


-20 . 0 


0.0 




J 


4 


-33.0 


0.0 




J 


5 


-100.0 


0.0 




J 


ZEROES 


REAL 


IMAGINARY 






1 


-5.0 


0.0 




J 


GAIN CONSTANT = 2640000 


.0 








TESTS FOR EXAMPLE 10 






RUN NO . 


OF POLES 


NO. OF ZEROES 


RESULT 


a 


5 


1 


See 


Table 


b 


3 


0 


See 


Table 


c 


2 


0 


See 


Table 



INPUT FUNCTION SIMILAR TO THE ONE SHOWN IN FIG. 5. 
INITIAL CONDITIONS EQUAL ZERO. 

SEE FIG. 6 COMPARING THE RESPONSES OF THEORETICAL 
AND IDENTIFIED SYSTEMS TO A UNIT STEP INPUT. 



Data 

Data 

Data 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




• 1 


-0.94505946 


o 

• 

o 


J 


2 


-10.807637 


o 

• 

o 


J 


3 


-32.034223 


o 

• 

o 


J 


4 


-62.693979 


o 

• 

o 


J 


5 


-120.97375 


o 

• 

o 


J 


ZEROES 


REAL 


IMAGINARY 




1 


-4.9340575 


o 

• 

o 


J 


GAIN CONSTANT = 2741280.0 






SYSTEM STATE 


VARIABLES (PHASE 


FORM) 




A VECTOR 








1 


248.538.7 






2 


2992971.5 






3 


403446.21 






4 


16013.272 






5 


227.45465 






B VECTOR 








5 


2741280.0 






C VECTOR 








1 


4.9340575 






2 


1.0000000 
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REMARKS: 

ASK FOR FIVE POLES ONE ZERO. 

RESULTS NOT VERY GOOD FOR THE HIGHER FREQUENCY POLES. 
SEE FIG. 6 COMPARING THE RESPONSES OF THEORETICAL AND 
IDENTIFIED SYSTEMS TO A UNIT STEP INPUT. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 



SYSTEM TRANSFER FUNCTION 
POLES REAL 

1 -2.6408186 

2 -24.059549 

3 -24.059549 



IMAGINARY 

0.0 

0.0 

0.0 



J 

J 

J 



ZEROES 



REAL 



IMAGINARY 



GAIN CONSTANT = 10410.082 



SYSTEM STATE VARIABLES (PHASE FORM) 
A VECTOR 

1 3331.3020 

2 1388.5395' 

3 50.759916 



B VECTOR 

3 10410.082 

C VECTOR 

1 1.0000000 



REMARKS : 

ASK FOR THREE POLES NO ZEROES. 

RESULTS CAN BE CONSIDERED A LOWER ORDER MATHEMATICAL 
MODEL. 

SEE FIG. 6 COMPARING THE RESPONSES OF THEORETICAL AND 
IDENTIFIED SYSTEMS TO A UNIT STEP INPUT. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (c) 



SYSTEM TRANSFER FUNCTION 

J 
J 

ZEROES REAL IMAGINARY 

GAIN CONSTANT = 182.89790 



POLES 

1 

2 



REAL 

-5.6608772 

-13.680660 



IMAGINARY 

0.0 

0.0 



SYSTEM STATE VARIABLES (PHASE FORM) 

A VECTOR 

1 77.444535 

2 19.341537 
B VECTOR 

2 182.89790 

C VECTOR 

1 1.0000000 

REMARKS : 

ASK FOR TWO POLES NO ZEROES. 

RESULTS CAN BE CONSIDERED A LOWER ORDER MATHEMATICAL 
MODEL . 

SEE FIG. 6 COMPARING RESPONSES OF THEORETICAL AND 
IDENTIFIED SYSTEMS TO A UNIT STEP -INPUT. 
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UNITS-OUTPUT 



25 



20 




0 12 3 4 



TIME-SECONDS 

UNIT STEP RESPONSE OF: 

(a) System to be identified 

(b) Fifth order identified model 

(c) Third order identified model 

(d) Second order identified model 



FIGURE 6 
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«• 

- 




<V 







EXAMPLE 11 





SYSTEM TO BE 


IDENTIFIED 




SYSTEM TRANSFER 


FUNCTION 






POLES 


REAL 


IMAGINARY 




1 


1 

H* 

• 

O 


o 

• 

o 


J 


2 


1 

O 

• 

O 


o 

• 

o 


J 


3 


-33.0 


o 

• 

o 


J 


4 


-100.0 


o 

• 

o 


J 


ZEROES 


REAL 


IMAGINARY 




1 


-5.0 


o 

• 

o 


J 



GAIN CONSTANT = 52500.0 

TESTS FOR EXAMPLE 11 

RUN NO. OF POLES NO. OF ZEROES RESULTS 

a 4 1 Good 

b 3 0 See Table of Data 

c 2 0 See Table of Data 

INPUT FUNCTION IS SIMILAR TO THE ONE SHOWN IN FIG. 5. 

INITIAL CONDITIONS EQUAL ZERO. 

SEE FIG. 7 & 8 COMPARING THE RESPONSES OF THEORETICAL 
AND IDENTIFYING SYSTEMS TO A STEP INPUT. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (a) 



SYSTEM TRANSFER FUNCTION 



POLES 


REAL 


IMAGINARY 




1 


-1.0182594 


o 

• 

o 


J 


2 


-10.853145 


o 

o 


J 


3 


-35.841486 


o 

• 

o 


J 


4 


-108.22539 


o 

• 

o 


J 


ZEROES 


REAL 


IMAGINARY 




1 


-5.0491588 


o 

# 

o 


J 



GAIN CONSTANT = 64711.781 



SYSTEM STATE VARIABLES (PHASE 
A VECTOR 



1 

2 

3 

4 

B VECTOR 
4 

C VECTOR 
1 
2 



42867.599 

47640.814 

5600.2860 

155.93828 

64711.781 

5.0491588 

1.0000000 



FORM) 



REMARKS: 



ASK FOR FOUR POLES ONE ZERO. 
RESULTS GOOD. 
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SEE FIG. 7 & 8 COMPARING THE RESPONSES OF THEORETICAL 



AND IDENTIFIED SYSTEMS TO A UNIT STEP INPUT. 



107 



















4 


















IDENTIFICATION OF UNKNOWN SYSTEM - Run (b) 



SYSTEM TRANSFER FUNCTION 
POLES REAL 

1 -31.063053 

2 -31.063053 

3 -2.9265019 



IMAGINARY 

-36.791868 

36.791868 

0.0 



J 

J 

J 



ZEROES 



REAL 



GAIN CONSTANT = 28337.477 



SYSTEM STATE VARIABLES (PHASE 
A VECTOR 



1 

2 

3 

B VECTOR 
3 

C VECTOR 
1 



6785.2550 

2500.3670 

65.052607 

28337.477 

1.0000000 



IMAGINARY 



FORM) 



REMARKS : 

ASK FOR THREE POLES NO ZEROES. 

RESULTS CAN BE CONSIDERED A LOWER ORDER MATHEMATICAL 
MODEL. 

SEE FIG. 7 & 8 COMPARING THE RESPONSES OF THEORETICAL 
AND IDENTIFIED SYSTEMS TO A UNIT STEP INPUT. 
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IDENTIFICATION OF UNKNOWN SYSTEM - Run (c) 



SYSTEM TRANSFER FUNCTION 
POLES REAL 

1 -3.4314173 

2 -36.703534 



IMAGINARY 

0.0 J 

0.0 J 



ZEROES 



REAL IMAGINARY 



GAIN CONSTANT = 537.30396 



SYSTEM STATE VARIABLES (PHASE FORM) 

A VECTOR 

1 125.94514 

2 40.134951 
B VECTOR 

2 537.30396 

C VECTOR 

1 1.0000000 

REMARKS : 

ASK FOR TWO POLES NO ZEROES. 

RESULTS CAN BE CONSIDERED A LOWER ORDER MATHEMATICAL 
MODEL . 

SEE FIGs. 7 & 8 COMPARING THE RESPONSES OF THEORETICAL 
AND IDENTIFIED SYSTEMS TO A UNIT STEP INPUT. 
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UNITS-OUTPUT 




TIME-SECONDS 
FIGURE 7 
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TIME-SECONDS 



FIGURE 8 
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IV. CONCLUSIONS 



The method of multiple integrations is a practical method 
for identifying lumped parameter, linear, time invariant 
systems. It works well in the absence of noise up to a fourth 
order system if the input function is carefully selected. 

For higher order systems the identification is in general not 
good, but an approximation of the system can be obtained. 

If desired, the method can also be used to obtain a reduced 
order mathematical model of the system, after the 
identification . 

The accuracy of an identification depends very strongly 
on the noise present. Small R.M.S. values of noise are 
enough to make the identification bad. Noise with a R.M.S. 
value of about 0.5 volts and below do not practically affect 
the identification, when inputs like the ones ued in this 
research are present, independently of where the noise is 
present - input, output or both -. Accuracy also depends 
very much on the type of input function used to drive the 
system. 

The computational requirements of the method are not 
excessive, however, an accurate A/D converter system is very 
desirable and the better this operation is done, the more 
accurate will the identification be, since the accuracy of 
an identification is usually comparable to the number of 
exact significant digits in the input-output data. 
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The limitation of the method to low order systems is 
primarily due to the algorithm used to solve the over-deter- 
mined set of linear equations. As better algorithms become 
available, it will be possible to identify higher order 
systems . 

The results obtained when reduction of the system order 
was desired were not very good at first sight, but when the 
steady state gain of the lower order models was adjusted, 
making it equal to the one of the non-reduced identified 
model, the approximation can be considered good. However, 
when a pair of dominant complex roots is present in the 
system a better reduced order model should be expected. 
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APPENDIX I 



A. OBTAINING THE INPUT - OUTPUT RECORD 

The identification technique described in Chapter II is 
based on knowledge of the input and output record and since 
the IBM 360 Digital Computer was used to perform the inte- 
grations of equation (5) , the record should be available in 
a digital form. Thus, and A/D (analog to digital) operation 
should be performed. This A/D operation was accomplished 
using the Hybrid Computer System (XDS 9300/Ci 5000) located 
in the Electrical Engineering Computer Laboratory. The 
output of this computer is a seven track magnetic tape in 
which the input and output functions are recorded in digital 
form. Since the XDS 9300/Ci 5000 system is used on a self- 
service basis, a brief discussion of what was done to obtain 
the record, is given here, so that the same procedure could 
be used in the future whenever an A/D conversion is needed - 
it should be pointed out here than the computer center staff 
was always predisposed to assist in the operation of the 
system. 

The analog part of the system, the Ci 5000, has two patch- 
boards which are the heart of the system. The logic patch- 
board was utilized to set the sampling rate, at which the 
analog signal was sampled - 10KHZ in this particular case. 

This sampling rate is the maximum that can be achieved 
with this computer, and can be changed (decreased) by 
altering a counter-like number display located in the 
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frontal panel of the computer - FPOO and FP04. Also in this 
patchboard the base frequency - 1KHZ - used to determine the 
moments at which the input function is changed step-wise, is 
set. The same counter-like number display can be used to 
change this frequency - FF10 and 14. In the digital patch- 
board, the wiring that controls the digital part of the A/D 
operation is also set. 

The analog patchboard was used to introduce the proper 
isolation between the different components of the different 
systems tested and to provide the desired amplification 
factor. This patchboard was also used in the generation of 
the input signal. This signal was obtained in the following 
way. An oscillator generates a sine wave in which amplitude 
and frequency can be adjusted by two potentiometers - P043 
and P417 -; the .sine-wave is amplified and passed through a 
zero order hold circuit which output is a step-wise wave form, 
that was used as the input function. The sine wave v/as also 
used to control the beginning of the sampling. The wiring of 
both patchboards was set by the computer center staff and 
those patchboards - numbered 21 and 21 - are reserved for 
the A/D operation. 

The input and output signals are fed from the analog 
part of the computer system to the digital part. The 
Digital part interfaced with the analog is a Xerox Data 
System - XDS9300. 

Two tapes drive units are available as input-output 
devices. One printer is used for program listings and 
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21. With the sense switches at the XDS9300 console off, 
depress CLEAR and CLEAR FLAGS buttons simultaneously. 

22. Mount the magnetic tape following the instructions 
available on the inside of the protective doors. 

23. Select density 556 bits per inch with the correspond- 
ing knob on the face of the tape driver. 

24. Read in the A/D program, described in Chapter III 

and listed in Chapter IV, following the instructions available 
on the console of the XDS9300. 

25. Enter, through the teletype, the parameters to be 
used for the A/D program. This program has seven options or 
instructions which are entered also through the teletype by 
a single digit number and depressing the carriage return 
button. If a further response is required the computer would 
indicate the format. The options are: 

1 = Desire to enter new parameters. 

2 = Start digitizing the analog input signals 
(normally actuated by manual switch operation on Ci 5000 
control panel; (DS1 up and key COMPUTE depressed) . 

3 = Write an END OF FILE on the magnetic tape. 

4 = Rewind the tape to the load point (beginning) . 

5 = Skip the following number of END OF FILE marks. 

6 = Print the following amount of digitized data. 

7 = Actuate the A/D program for the next single 
block of data encountered and feed the signal to the strip 
chart recorder (v/hich must be running prior to this option 
selection connected to the analog patchboard on the Ci 5000. 
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The rest of the counters should be set to zero. 

7. Depress the master clock control, RUN. 

8. Set DFOO (Delay flip-flop) to .1 ms. and adjust it, 
until its corresponding light on the front panel of the 



Ci 5000 


comes on. 


9. 


Set DF06 in .1 ms., DF07 in 1 sec. and DF10 in 1 ms. 


10. 


Set all the DS (digital switches) to the center 



position. 



11. 


Turn on the TRC . 


12. 


Set switches XI to the left and X2 to the right. 


13. 


Select amplifier 042 in the Ci 5000 console, adjust- 



ing P417, DF06 and F10, F14 to obtain the desired input 
function which can be observed in the TRC when the DS1 
switch is set in the up position and the button COMPUTE 
is depressed in the Ci 5000 console. 

14. Set the R.T . clock switch on XDS9300 to EXIT. 

15. If the XDS9300 is not energized, depress simultaneously 
RESET and POWER buttons. 



• 

kO 

i— 1 


If the teletype is not ON, energize it by turning the 


knurled 


knob to ON. 


17. 


If hte Line Printer is not ON, depress the POWER 


switch • 


- indicator on the face of the printer. About one 



minute later the lower half of the Power button/lights up 
indicating the printer is energized. 



• 

00 
1 — 1 


Depress READY SWITCH INDICATOR. 


19. 


Depress the POWER ON button on the card reader. 


20. 


Depress IDLE button on the console of the XDS9300. 
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outputs. A card reader is associated with the digital com- 
puter as the avenue of input for compiling programs. A 
teletype unit is used for parameter . inputs and instructions 
to the digital computer program. A control console displays 
information of the computer and has the corresponding control 
switches. A recorder is also available with which the 
records of the input and output, after the A/D operation is 
finished, can be recorded in a strip chart. 

After the preceding general description of the system, 
the following procedure can be followed to take a record. 

1. Set the analog and digital patchboards in the 
Ci 5000. 

2. Connect the output of the amplifier 042 to the 
input of the system. 

3. Connect the output of the system to the input of 
the amplifier 036. 

4. If the Ci 5000 is not energized, depress on switch 
indicator at lower left corner of display panel of Ci 5000. 
Depress KEY BOARD and POTSET buttons on the display panel. 

5. Enter the maximum value of the sine wave through 
the POTSET, POT and SERVO keys. 

6. Set the counters FF ' s in the following positions: 

FF00 in 0, FF04 in 9, for 10KHZ sampling rate. 

FF10 in 2, F14 in 4, for 40HZ rate of changing, 

step-wise, the input function. 

FF30 in 0, if only one record per run is desired; 
increasing FF30 the number of records/run can be adjusted. 
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The parameters to be entered in point 25 are: 

NSAMP = Total number of samples to be taken from 
each channel - 4002 in this case. 

NCHAN = Number of channels to be digitized - 2 in 
this case. 

NREC = Number of records to be digitized per run; 
this number should be equal to the one set in the counter FF30 
(Ci 5000) plus 1; in this particular problem NREC = 1. 

ITAPE = The identifying number of the tape drive in 
use. This number was set when the tape was mounted. 

NDEL = Controls the delay between the digital and 
analog signal in option 7. 

After entering the parameters the * character should be 
typed on the teletype followed by RETURN. At this point 
everything is ready to take a record; by typing, 2 - RETURN, 
on the teletype, and setting switch DS1 - UP and pressing 
COMPUTE on the Ci 5000. After the record was taken, write 
- 3 - (END OF FILE). If it is desired, the record can be 
recorded in a strip chart recorder. 

B. SEVEN TO NINE TRACK TAPE CONVERSION 

In the magnetic tape generated in point A, the input- 
output record is written in a seven track tape, but since 
this tape is going to be used with the IBM 360 Digital 
Computer and this computer uses NINE TRACK TAPES, a 
conversion of the seven track tape into a nine track tape 
is required. 
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This operation is done by the IBM 360, using the program 
CONVERT, which is described in the Technical Note No. 0211-08 
second edition - PROCEDURES FOR CONVERTING 7 -TRACK MAGNETIC 
TAPES TO 9 -TRACK MAGNETIC TAPES by Sharon D. Raney [21]. 

This description is available at the Computer Center in 
Ingersol Hall. A listing of the CONVERT program is given 
in Chapter IV. 
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HIS PROGRAM is 78 BE USE-) WITH THE C I 5000- SDS9300 HYBRID COMPUTE 
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THIS PROGRAM CONVERTS A SEVEN TRACK TAPE TO A NINE TRACK TAPE 
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THIS PROGRAM SOLVES A SET OF LINEAR DIFFERENTIAL EQUETIONS 
THE SOLUTION APPEARS AS A DRAW 
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